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Abstract 

Satellite observation of specific Earth sites or regions usually 
requires the use of compatible (resonant) orbits. If the 
observing sensor is optimized for a given altitude then the 
orbit must also be circular. This study provides analytical 
derivation of satellite ground tracks of any compatible 
circular orbits and for a full repetition cycle. The purpose is 
to offer equations whose solutions are all the points 
coordinates in which the ground track crosses itself, here 
called -points. The geographical locations of these points 
along with the approximate expressions of the associated 
dwell times allow developing optimal orbit design for 
civilian or military satellite observation missions. 
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Nomenclature 

a= semi major axis [km] 
e= eccentricity [-] 
p=semi latus rectum [km] 
;'=inclination [rad] 

/2=right ascension of the ascending node [rad] 

<y=anomaly of perigee [rad] 

(p= true anomaly [rad] 

n=mean motion [rad/s] 

/l=latitude [rad] 

P= longitude [rad] 

X P =orbital nodal period [s] 

T.i=Greenwich nodal period [s] 

/2=Earth's second dynamic form factor [-] 

R©=Earth's radius [km] 

oa® =Earth sidereal rotation rate [rad/s] 

/^Earth's gravitational constant [km 3 /s 2 ] 

Cx=cos x 

Sx=sin x 

Introduction 

Observation by satellites on Earth's sites or regions is 



of interest for both civilian and military purposes. 
Space missions to monitor the values of physical 
quantities in specific geographical areas are of 
great interest for the scientific community while 
reconnaissance missions to perform persistent 
observation on sites, borders, or regions, are of high 
interest for intelligence and the defense community. 
Another important area of interest is satellite 
communications with specific Earth locations by single 
satellites or by satellite constellations. 

Observation on Earth's sites or regions was also the 
focus of «Two-Way Orbits», from O. Abdelkhalic e D. 
Mortari, where the concept of " two-way s" orbits was 
introduced. The repeated ground track path of these 
orbits is a closed-loop trajectory that intersects itself, in 
some points, with tangent intersections. This allows 
designing orbits (and constellations) where the 
spacecraft passes over these tangent intersections once 
in a prograde and once in a retrograde mode and at 
different altitudes. The design of these orbits has been 
done by numerical optimization, specifically using 
Genetic Algorithms. Conversely, the focus of this 
paper is to find the analytic expressions of the 
intersection points of repeating ground-track 
trajectories (called X-points) as produced by circular 
orbits. Motivation is driven by the fact that the use of 
elliptical orbits is in contrast with the design of 
observing sensors, which is optimal for specific fixed 
altitude. The design of orbits that optimally and 
repeatedly observe Earth's sites implies the use of 
compatible orbits whose definition is summarized in 
the next subsection. 

Compatible Orbits 

To obtain optimal coverage of specific Earth's sites 
or regions, compatible orbits must be adopted. 
According to Carter, these are a set of orbits for which 
the ratio between the orbit nodal period (time between 
two subsequent equator crossings in ascending mode, 
T P ) and the nodal period of Greenwich (elapsed time for 
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Greenwich meridian between two subsequent 
crossings of the orbit ascending node line, Td) is 
rational. Carter defined these orbits as 

N p T p = N p ^r = N d T d = N d ^- (2.1) 
n + co co^-Ll 

where N P and Nd are two integers representing the 
number of nodal periods completed by the satellite 
and the number of times that the Greenwich meridian 
crosses the orbit ascending node line in the repetition 
time (cycle), T rep =N P T P = NdTd, respectively. The main 
gravitational perturbation has been taken into account 
in Equation (2.1) due to the Earth oblateness, known as 
the ]i effect. The secular and persistent ]i effect 
modifies the mean motion according to 

R&, 

(2.2) 



1 + -/. 



(2-3sin 2 oVl-e 2 

4 { p 

where no =^JJtJa? , the unperturbed mean motion, 
linearly changes the right ascension of the ascending 
node, 

as well as the argument of perigee, 
4 \ P 

while the remaining parameters (a, e and i) remain 
unchanged 

da de di 
— »0, — «0,and — «0. 
dt dt dt 

In particular, for circular orbits (e = 0) the sum, n + co, 
assumes the value 

= n + cb = n [1 + £(2 - 3 sin 2 /)] [1 + £(5 cos 2 i - 1)] 

where 



ncos (, 



«(5 cos 2 /-!), 



(2.3) 



(2.4) 



4 : 



However, by disregarding the f term (for R = 8,000 
km, ^= 5 x 10 4 ), previous equation can be written as 

n c = n + d> «n [1-2^(1- 4 cos 2 0], (2.5) 

which is an approximated expression for the modified 
mean motion of circular orbits. 

The orbit nodal period as defined by Eq. (2.1) is correct 
for circular orbits only. For elliptical orbit the nodal 
period is also a function of the initial value of the 
argument of perigee, coo. However, the variation of the 
nodal period with a>o is very small (order of few 
seconds for LEO orbits) and therefore, Carter's 
definition, T P =2n/(n + oj) can be considered a good 
expression of the nodal period for small eccentric 
orbits. The rotation of the apsidal line affecting 
elliptical orbits complicates the definition of 
compatible orbits as the apsidal line rotation period 



must be synchronized with the orbit nodal period. 
For assigned values of Nd, N P , e, and i, Eq. (2.1) 
becomes a nonlinear equation in terms of the semi- 
major axis only. This nonlinear equation must be 
solved to obtain the compatible semi-major axis. This 
is done by rewriting Eq. (2.2), Eq. (2.3), and Eq. (2.4), in 
the following compact mode 



CO = c„. 



(2.6) 



Where 



l-e 1 



(2-3sin 2 /)Vl-e 2 



c n = 



c„, = 



2 \\ 



Rc. 



cos; 



4 2 {l-e 2 



(5cos 2 /-l) 



Substituting Eq (2.6) with Eq (2.1) and setting c n = 
N p c n + N d c M , the following function, 



ffa) = N(m m -Cl)- 



a \ a 



1 + - 



: 0, (2.7) 



is obtained. Equation (2.7) can be solved using 
Newton-Raphson iterative method, a k+1 = a k — jr^^, 

where, by setting c = c m N d + c n , the derivative has 
the expression 

f(a) = ^3N d+ 7^ + U i - 
2a\ a 

and where the starting point ao can be chosen as the 
semi-major axis of a compatible unperturbed orbit 



4 

a j 



N„co m 



N„ 



N p co e 



(2.8) 



The report «Goddard Trajectory Determination 
Subsystem Mathematical Specifications^ by Wagner 
presents an alternative iterative method to calculate 
the mean semi-major axis of a compatible orbit. The 
algorithm begins with the evaluation, as a starting 
point, of the unperturbed semi-major axis ao, as in Eq. 
(2.8). Then, a is evaluated through the following 
iterative process 



where 



1 + -/J — 

4 a b 



(2-3sin 2 



2N „ 
1 H -cos, i 



-5 cos i 
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In order to obtain the equation of the ground track, the 
next subsection provides the trigonometric 
relationships relating the angles used to identify the 
position of a point on the Earth surface (geographical 
latitude, A, and longitude, P), and the angular 
parameters used to define the position of the satellite 
in the inertial frame (inclination, i, and the anomaly 
from the node line, <p). These relationships are 
introduced and the ground track parametric equation 
is then derived. 

Geometrical Relationships 

Let us consider the ground track of a satellite in 
circular orbit with inclination Without losing 
generality, it can be assumed that the right ascension 
of ascending node is aligned with the first point of 
Aries (vernal axis), Q = 0, and the Greenwich meridian 
crosses the ascending node, «g = 0. 




FIGURE 1 SATELLITE OBSERVATION GEOMETRY (ASCENDING 
MODE) 

Under these assumptions, it can be seen from Fig. 1 
that the angles involved for the ascending point 
(subscript "a") are 8 a , indicating the right ascension 
from the Y axis, A, the geographical latitude of the sub- 
satellite point, and q> = nd, indicating the angle in the 
orbit plane from the sub-satellite point to the 
ascending node line. Among these angles, the 
following relationships hold 



cos i = sin % a cos A 



tan 8 . = cos i tan < 



tan i tan % a = 1 / cos </> a 



cos x a = sm i cos S a and <^ sin 8 a = sin % a sin $ a (2.9) 



sin i sin <f> a = sin A 



All these relationships are just spherical trigonometry 
identities applied to the spherical triangle which is 



highlighted in Fig. 1. For the descending point 
(subscript "d"), which happen at the same latitude, the 
angles satisfy the following relationships 

<l>a= n -<t>a> S d =n-S a , Xd =n-x a - 
Ground Track Equation 

The parametric equations for the ground track can be 
obtained considering the unit vector identifying the sub- 
satellite point in the Earth-Centered Earth-Fixed (ECEF) 
reference frame, s E = {cos A cos/?, cos A sin/?, sin A} r , 
where A and /? indicate the sub-satellite point latitude 
and longitude, respectively. The position of the 
satellite in the Orbit Inertial frame is given by the unit 
vector s s = {cos n c t , sin n c t,0} T . For circular orbit and 
disregarding the effect of Earth oblateness, the 
transformation matrix between these two reference 
frames is R 3 (oj@t — fl)R 1 (t), where Ri and R3 are the 
matrices performing rigid rotations about the first and 
third coordinate axis, respectively. Therefore, we can 
write 



c c 




~c a 
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where a = — Q)t which leads to the following 
system of equations 



c c 

A p 

CxSp 



■c a c i c t + ss t 



1 - S { S nt 

that allows us to obtain the equations of the ground 
track 




Figure 2 shows, by setting Np = 24, Nd = 7, and i = 55°, 
an example of the ground track trajectory for the 

Ik Ik 

repetition time T = A 7 " — = N d ■ in which two 

"elementary patterns," associated with the north and 
south hemisphere, respectively coloured in red and 
green are highlighted. 

In addition to the obvious and expected symmetry 
with respect to the axes, several regularities are visible 
in the pattern of the ground track trajectory. In 
particular, for a given value of latitude, there are 
exactly N P intersections, equally spaced in longitude. 
Similarly, for an assigned value of longitude, several 
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intersections at different latitudes occur. The two 
parameters ruling the intersection pattern are the 
radius (circular orbits) and the inclination. In the next 
subsection the conditions capable to evaluate the 
ground track intersections are given. 

N = 24, N = 7, i = 55 deg. Number of intersections = 384 




-150 -100 -50 50 100 150 

FIGURE 2 EXAMPLE OF SATELLITE GROUND TRACK 
INTERSECTIONS DURING A REPETITION PERIOD 

Solving the Compatibility Condition 

Let us consider a point where the ground track 
trajectory crosses a given latitude for two times, one 
when it ascends toward north and another time when 
it descends toward south. Let AT be the time 
difference between the time at which the satellite 
crosses that latitude in ascending mode and the time 
when it crosses again in descending mode. This time 
interval can be expressed in two different ways, using 
either Earth rotations 



2na S A — 8„ 



2na n - 28, 

- + a eU : 



AT 

co 9 - Q co m - Q a> m - Q (o m - Q. 
or orbit revolutions 

2np <f> d -0 a _2np n-2<j> a 



(3.1) 



AT = - 



- + - 



- + - 



^ eD - (3.2) 



Both these equations account for an integer number of 
complete rotations or revolutions plus the exceeding 
angles &+& or (p&+(p*. Employing Equating Eq. (3.1) 
and Eq. (3.2) as well as the compatibility equation for 
circular orbits, N p (a> @ — fi) = N d n c , we obtain 

— ^(2xa + x-2S a ) = —(2xfi + x-2fi ) 
(O m -n n c 

which can be written as 

. 7 In In n-2(j) 

aN n + N J> - = /? — + ^. 

N.n N.n n n 

a c d c c c 

Using the third of Eq. (2.9) and considering that cp^nd, 
we obtain 

^■-2tan"'[cos! tan(n.f)] », n — 2nJ. T 
aN p -pN d + -^-N,, : L -N„ 



2/r 



2tt 



(3.3) 



The first two terms in Eq. (3.3) always add up to an 
integer. The last two instead can take any real value 
and can be indicated by a function of time/ft). When/ 
is an integer, then Eq. (3.3) becomes a Diophantine 
equation whose solutions identify intersections of the 
ground track trajectory. The analysis of the /function 
is the main focus of this study. Because of the 
symmetry, the analysis can be limited to the n c t £ 
[o, |] range. 

Analysis of Function f 

The purpose of our analysis is to find the times t at 
which the function / assumes integer values. The 
function is continuous in nd and, because of the 
symmetry, the search can be limited to the n c t £ [o, |j 
range. Therefore, finding minimum and maximum 
values of / in the searching interval, called [fmin,/max], 
also provides the complete set of integers for / and 
associated times where the ground track intersects 
itself, i-e.[\f min ],[f max \] El r M and |»J indicate the 
"floor" and the "ceil" functions, respectively. 

These extremal values may occur at the range limits 

(nd= and nd = n/2) or inside the interval, in which 

case the first derivative /' can be used to find the 

locations of minumum/maximum (we will show later 

that /can have at most one internal stationary point). 

Since f(0) = (N v - Nd)/2, there will be a set of 

intersections at the equator (nd = 0) whenever JVp and 

Nd are both odd (they cannot be both even as they 

must be coprime), so that their difference is even. As 

nd approaches n/2, function /assumes different values 

depending on the sign of cosf because 

f+;r/2, if cos / > 
lim tan (cos i tan x) = < 
*_►£ \-nl2, if cos i < 

2 k 

and, in particular, f(n/2)= for prograde orbits and 
f(n/2) = N P for retrograde orbits. 



The expression of the first derivative,/', is 

df l 



N., 



N cos; 

v 



B(n c t) 

while the second derivative,/", is 

a 2 / _ n i 



cos («.f) + cos /sin (n.t) 



(3.4) 



d(n c tf 



cos (« t .r) + cos /sin (nj) 



cos/sin(2n ( ./)- (3-5) 



aN-/3N d +f{t) = Q. 



Note that /" is always negative for prograde orbits 
(cos i > 0) and positive for retrograde orbits (cos / < 0), 
which implies that/' is monotone in all cases, proving, 
as previously anticipated, that / can only have one 
internal stationary point (minimum or maximum). The 
analysis so far clearly shows how the inclination / 
plays a primary role in determining the behavior of /. 
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Its importance will be highlighted further in the 
discussion regarding the zeroes of/'. 

By setting /"= 0, we obtain 

sm 2 (n c t)cos 2 i-TCOsi + cos 2 (n c t) = 0, (3.6) 
where r = N p I N d , which can be reformulated in two 
different forms 

l — T cos i 



sin (nt) = - 



sin 2 i 



(a) 



cos J 



-■sjr 2 - sin 2 (2nt) 



(b) 



(3.7) 

2 sin (n c t) 

where the positive solution of Eq. (3.7b) is discarded 
as it gives values for cos; outside the admissible range 
Equation (3.7a) enable to derive the limit 
values of inclination for which the derivative can 
vanish 



sin (n c t) = -> z min = cos" 



sin 2 («.?) = ! -> / max = 



(3.8) 



Therefore, there is no maximum within the [/(0),/(jt/2)] 
range for orbits at i < cos -1 (1/x). Another way to obtain 
this result is by means of the note that i < cos -1 (1/x) is 
equivalent to (Np/Nd) cos;> 1. By rewriting Eq. (3.6) in 
the form (sin 2 i cos 2 (nd) + cos 2 i) -(N P /Nd)cos i, and 
since the terms in brackets cannot be greater than 1, it 
is concluded that the derivative assumes negative 
values at both extremes. Together with monotonicity, 
this proves the result. 

In the following, we indicate cos -1 (1/x) = ii, and this 
inclination is called the first characteristic inclination. 
When i = ii, the maximum is located exactly at nd = 
with/' = 0. Orbits with inclinations i < n "experience the 
maximum with /' < 0. Assuming the spacecraft at the 
ascending node at the initial time then / is a 
representation of the first quarter of a ground track. 
Consequently, an orbit inclined at the first 
characteristic inclination, ii, crosses the equator along 
the local meridian. Whenever N P - Nd is even, 
subsequent ground tracks are also tangent to each 
other at the equator. 

Alternatively, Eq. (3.7b) can be substituted into Eq. (3.5) 
to just obtain a function of nd 

T-^jr 2 -sin 2 (2n f f) | 

2sin(2n ' (3-9) 



N-N 



Equation (3.9) provides solutions (associated with 
/* E Z) of {nd} for which the ground track is tangent to 
itself at those phase values. Once substituted in Eq. 



3.7b), these {nd} values identify the characteristic 
inclinations that give the tangency condition (the 
special case for nd= can be regarded as a limit). These 
constitute a set of characteristic inclinations 
determined by the values of N P and Nd, andtherefore, 
they are specific for every orbit. Equation (3.7b) also 
implies that the characteristic inclinations approach 
90° monotonically as nd increases. Whether a certain 
compatible orbit has one or more characteristic 
inclinations depends on the value of r. Empirical 
analysis shows that there will always be at least 2 
characteristic inclinations as long as t > 5/3. 
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FIGURE 3 CRITICAL INCLINATIONS FOR A COMPATIBLE 
ORBIT. 

As previously discussed, subsequent ground tracks of 
an orbit inclined at ii with N P both and Nd odd will be 
tangent at the equator. If the inclination is equal to ii, 
then subsequent ground tracks will be tangent at the 
associated phase value. In general, an orbit at the iV 
characteristic inclination (with k>2), will have its first- 
passage ground track tangent to the one corresponding 
to the k-1 subsequent passage, respectively. 

For retrograde orbits (i > 90°) the first derivative is 
positive at both extremes of the interval. Therefore, no 
stationary points are possible for / which takes strictly 
monotonically increasing values in the interval [(N P - 
Nd)/2, N P ], and the ground track of a retrograde orbit 
will never be tangent to itself at any point. 

In conclusion, / assumes integer values in the interval 
[fa, fi>], according to the following scheme (for the 
unique case of N P = 3 and Nd = 2, ii is not defined and 
fi*has to be used in the formulas instead): 



< i < L 



• L <i<nll 



i>xl2 



0, 



/ e [o,[ 



Z(V)J] 



N-N, 



N. 
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For clarity sake, the results of our analysis are 
summarized here: 



1. 



Function / assumes integer values in the 
interval [[/ m ;„l, [f ma x J]- F° r prograde orbits, /mm 
= 0, while for retrograde orbits /mm = (Np - Nd)/2. 

As long as i < if, function / is monotonically 
decreasing. The greatest integer is equal to 
[(.N p -N d )/2\. 

If Np - Nd is even at the first characteristic 
inclination, the ground track goes along the 
local meridian (perpendicular to the equator). 
Subsequent ground tracks are also tangent 
there, and for greater inclinations the ground 
track shows a "reversed" intersection pattern 
compared to the typical one of lower 
inclinations, but which keeps the same 
properties of symmetry and pattern repetition. 

If Np - Na is odd, then subsequent ground 
tracks will be tangent to each other for the first 
time at the second characteristic inclination, at 
the phase (nd)*, 



(nt) =■ 



1 



sin i. 



N„ 



-COS J, 



(3.10) 



that always falls inside the searching interval, 
and the reversed pattern appears for 
inclinations greater than tt. 

5. As nd increases, function / assumes one or 
more integer values. For each of these values 
there will be 2 sets of N v intersections 
symmetric with respect to the equator, one 
located along the orbits at phase (nd)", and the 
other at phase n - (nd)". If the orbit is also 
inclined at one of the characteristic inclinations 
associated with the (nd)" values then at (nd)" 
the ground track intersects itself with vertical 
tangent. 

6. For retrograde orbits, function / has no 
stationary points within the interval. The 
maximum occurs at the right extreme of the 
interval, where the value is Np. 

Dwell Time Estimation 

A first order approximation of the dwell (observation) 
time of the satellite above a certain intersection point 
can be obtained. Given the orbit inclination i, the 
radius r, and setting simplicity Q = 0, the velocity 
vector of the satellite in the ECI frame can be easily 
written as 



v„, = 



— sin(n c t) 
cos(n c t) V 




In the ECI frame, the Earth rotational velocity at the 
given intersection point is 



-sm(a> B t) 
=1 ra e x(R s ?,.)!<{ cos(ty e f) 




= R s a> m cosA 



-sin(<» e f) 
cos(a) e f) 




Therefore, the relative velocity between satellite and 
intersection point can be obtained by simply adding 
the two vectors 

-sin(n c f) 
cosi'cos(n,f) 
sin / sin(n c f ) 



■ R s a m cos X - 



-sin(» e r) 

COS(ffi> ffi f) 





The magnitude of this vector is 
v 2 =vL+vl t [l + cos(2n.t)] + 
2v site v sat [sin(ftx B 0sin(n t) + cos icos(co e t)cos(n t)] 



where Vsite=R m co m cosA and v sat = Jn/r. As 



first 



approximation, this relative velocity can be assumed 
constant during the whole observation time. Therefore, 
given the nominal altitude of the satellite, h = r - R m , 
the observing sensor field-of-view angle, a, and 
ignoring the curvature of the Earth, the distance 
traveled by the satellite while in view of the target can 
be linearly estimated and consequently, the 
observation dwell time can be computed, as shown in 
Fig. 4. 




FIGURE 4 SATELLITE PASSAGE ABOVE AN INTERSECTION 
POINT 

In fact, the geometry allows us to write 

, „ , . 2d 2h sin a 
dUhsma -> t AmSl = — = . 



Numerical Examples 

In the following section we illustrate different possible 
cases associated with the various behaviors of the 
parameter / as i increases. For each case, the 
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intersection pattern both in a 3-D and 2-D view is 
shown, plus the plot of the parameter / and the 
passage above a certain intersection point with the 
associated dwell time. 

Case 1 - i<ii* 

This is the case relative to a prograde orbit whose 
inclination is lower than the characteristic value. The 
behavior does not change if the sum N P + Nd is even or 
odd. The repetition pattern for a resonant orbit with 
i < ii and ground track are shown in Fig. 5 and Fig. 6, 
respectively. Figure 7 shows the behavior of the 
parameter / and Fig. 8 shows the passage of the 
satellite above an intersecting X-point. 




-ISO -120 -SO 60 120 ISO 



FIGURE 6 2-D VIEW 



T 




1& 20 3O 40 5OW70 90M 



FIGURE 7 PARAMETER f FOR A RESONANT ORBIT WITH i< ii* 



Dwell lime above point {-65, 65): 441 s 




46 50 60 70 BO SO 



FIGURE 8 DWELL TIME FOR A RESONANT ORBIT WITH i< ii* 

Case 2 - Np+Nd even, i sii* 

This is the case of an orbit with even N P + Nd, and 
inclination close to the characteristic value h*. The 
repetition pattern for a resonant orbit with even Np + 
Nd at characteristic inclination is shown in Fig. 9 and 
Fig. 10, where it is possible to observe how successive 
tracks are tangent at the equator. Figure 1 1 shows how 
the parameter / has a maximum at nd = 0. The dwell 
time is shown in Fig. 12. 
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FIGURE 10 2-D VIEW 
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FIGURE 11 PARAMETER f FOR A RESONANT ORBIT WITH N P + 
Nd EVEN, AT CHARACTERISTIC INCLINATION 

&*ell lime Above BOtfri IflO, 57): 163 a 




20 36 40 50 60 70 BO M 



FIGURE 12 DWELL TIME FOR A RESONANT ORBIT WITH N P + 
M/EVEN, AT CHARACTERISTIC INCLINATION 

Case 3 - Np+Nd even, i=ii 

This is the characteristic case for all the orbits with odd 
Np + Nd. Repetition pattern for a resonant orbit with 
odd Np + Nd at characteristic inclination, is shown in 
Fig. 13 and Fig. 14, where the tracks are tangent not at 
the equator, but at two specific latitude values 
symmetric with respect to the equator itself. It is 
possible to notice from Fig. 15 how the parameter has 
its maximum equal to an integer. 




FIGURE 13 3-D VIEW 

N ■ 13, N ■ 4. i ■ as dea Number or interaectwrs ■ 1M 
p s 




FIGURE 14 2-D VIEW 
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FIGURE 15 PARAMETER / FOR A RESONANT ORBIT WITH Np + 
Nd ODD, AT CHARACTERISTIC INCLINATION 
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FIGURE 9 DWELL TIME FOR A RESONANT ORBIT WITH Np + Ni 
ODD, AT CHARACTERISTIC INCLINATION 

Case 4 - i>i" 




FIGURE 17 3-D VIEW 
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FIGURE 10 PARAMETER/FOR A RESONANT ORBIT AT 



INCLINATION GREATER THAN CHARACTERISTIC 
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FIGURE 11 DWELL TIME FOR A RESONANT ORBIT AT 
INCLINATION GREATER THAN CHARACTERISTIC 

When the inclination goes beyond the characteristic 
value, the satellite will move opposite to the Earth 
rotation, which causes a change in the pattern of 
intersections. The repetition pattern for a resonant 
orbit at supercritical inclination is shown in Fig. 17 and 
Fig. 18. The maximum of / increases further and 
reaches higher integer values, thus new intersections 
appear (Fig. 19). 

Case 5 - i>90° 

When the orbit is fully retrograde, the parameter /will 
become monotonically increasing (Fig. 23). The 
repetition pattern for a resonant retrograde orbit with 
different type of intersection pattern is shown in Fig. 
21 and Fig. 22. 




FIGURE 21 3-D VIEW 
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FIGURE 22 2-D VIEW 
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FIGURE 23 PARAMETER F FOR A RESONANT RETROGRADE 
ORBIT 
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FIGURE 24 DWELL TIME FOR A RESONANT RETROGRADE 
ORBIT 

Conclusions 

This study contains analytical equations developed for 
optimization of satellites and constellations for Earth 
observation using circular and compatible orbits. 
Based on intrinsic properties of a repetitive ground 
track, which presents many intersection points during 
a repetition period, the relations between the orbit 
parameters and the coordinates of the X-points, are 
revealed. Therefore, the orbit can be designed to take 
the biggest possible advantage from the innate 
symmetry of the intersection pattern. For instance, if a 
mission requires the satellite to have a X-points at 
certain specific coordinates, the other X-points can be 



used to enhance the scope of the mission, if they 
happen above other interesting sites. Conversely, if 
repeated coverage of a certain zone is required, 
consideration on the dimension of the field of view of 
the satellite is possible to design the orbit so that the 
zone is surrounded by the highest possible number of 
X-points which have coverage 2 to 4 times more 
frequent than that for usual orbits. 
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